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ABSTRACT 

We fit axisymmetric 3-integral dynamical models to NGC 3379 using the line-of-sight velocity distribution 
obtained from HST/FOS spectra of the galaxy center and ground-based long-slit spectroscopy along four position 
angles, with the light distribution constrained by WFPC2 and ground-based images. We have fitted models with 
inclinations from 29° (intrinsic galaxy type E5) to 90° (intrinsic El) and black hole masses from to 1O 9 M . 

The best-fit black hole masses range from 6 x 10 7 -2 x 10 8 M Q , depending on inclination. The preferred incli- 
nation is 90° (edge-on); however, the constraints on allowed inclination are not very strong due to our assumption 
of constant M/Ly. The velocity ellipsoid of the best model is not consistent with either isotropy or a two-integral 
distribution function. Along the major axis, the velocity ellipsoid becomes tangential at the innermost bin, ra- 
dial in the mid-range radii, and tangential again at the outermost bins. The rotation rises quickly at small radii 
due to the presence of the black hole. For the acceptable models, the radial to tangential ((og +ai)/2) disper- 
sion in the mid-range radii ranges from 1.1 < a r ja, < 1.7, with the smaller black holes requiring larger radial 
anisotropy. Compared with these 3-integral models, 2-integral isotropic models overestimate the black hole mass 
since they cannot provide adequate radial motion. However, the models presented in this paper still contain restric- 
tive assumptions — namely assumptions of constant M/Ly and spheroidal symmetry — requiring yet more models 
to study black hole properties in complete generality. 
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1. INTRODUCTION 

The dynamical state and mass distribution in the central re- 
gions of elliptical galaxies provide clues to the formation and 
evolutionary histories of these galaxies (see Merritt 1999 for a 
review). Consequently, the existence of central black holes has 
been the target of intense scrutiny (see Kormendy & Richstone 
1995 for a review). In this paper, we study the elliptical galaxy 
NGC 3379 using both Hubble Space Telescope and ground- 
based photometric and kinematic data in order to understand 
its central dynamics. 

At ground-based resolution NGC 3379 is a prototypical ellip- 
tical (El) galaxy with M v = -20.6 (Faber et al. 1997) at a dis- 
tance of 10.4 Mpc (Ajhar et al. 1997). Lauer et al. (1995) clas- 
sify it as a "core" galaxy — a galaxy that has a break in the sur- 
face brightness profile but still maintains a rising profile into the 
smallest measured radius. Faber et al. (1997) hypothesize that 
core galaxies are associated with massive black holes. How- 
ever, few core galaxies have strong black hole detections: M87 
(Harms et al. 1994), M84 (Bower et al. 1998), and N4261 (Fer- 
rarese et al. 1996). These detections result from gas dynam- 
ics, not stellar dynamics. Stellar dynamical evidence for central 
black holes in core galaxies is difficult to obtain since (1) core 
galaxies have low surface brightnesses, making kinematic ob- 
servations difficult, (2) the lack of rotation complicates the dy- 
namical modeling since velocity anisotropics potentially gov- 
ern the dynamical support, and (3) stars near the center travel 
out to the core radius so their contribution to the central velocity 
profile is de-weighted (Kormendy 1992). NGC 3379 was cho- 
sen for this study because of its relatively high central surface 
brightness among core galaxies, and because of its otherwise 
normal morphological and dynamical structure. 

The proximity, brightness, and normalcy of NGC 3379 
make it the object of many ground-based photometric (de Vau- 
couleurs & Capaccioli 1979, Lauer 1985, Capaccioli et al. 
1987, Capaccioli et al. 1991) and kinematic (Kormendy 1997, 
Kormendy 1985, van der Marel et al. 1990, Bender et al. 1994, 
Statler 1994, Statler & Smecker-Hane 1999) studies. However, 
previous models for NGC 3379 have yielded inconsistent re- 
sults. Assuming that the departure of the light profile from 
the R 1 ' 4 law signifies the existence of a stellar disk, Capacci- 
oli et al. (1991) find that the most likely model is nearly face-on 
with an inclination of 31° (intrinsic E5 galaxy). Their conclu- 
sion was that NGC 3379 and NGC 3115 have the same intrinsic 
shape but are seen from different viewing angles. In their model 
for NGC 3379 the central region is nearly oblate, but at about 
r = 0.5r e (r e = 35") the model becomes triaxial. No dynami- 
cal information was used in this model. A major uncertainty 
is that slight departures from a R l l A law, on which their model 
is based, may be due to tidal effects (Kormendy 1977) instead 
of a disk. Furthermore, there is no reason to expect galaxies to 
follow a R l l A law exactly. 

Kinematic information is required to determine unambigu- 
ously the inclination of a spheroidal system. Van der Marel 
et al. (1990), using the ground-based velocity data of Davies 
& Birkinshaw (1988) and Franx et al. (1989), conclude that 
the inclination of NGC 3379 is 60°; however the two datasets 
are not consistent and give slightly different velocity anisotropy 
parameters. Van der Marel et al. used 2-integral flattened mod- 
els with a parametrized form for the anisotropy. No goodness- 
of-fit criteria were given, and thus we are not able to compare 
their results with those of others. Statler (1994), using Bayesian 
statistics on the same datasets, prefers models that have incli- 



nation greater than 45°. Statler rejects Capaccioli's conclusion 
that the galaxy is flattened and triaxial at the 98% confidence 
level, but his results could accommodate those of van der Marel 
et al. Statler uses only the mean velocity measured along four 
position angles; the radial variation in the velocity or the dis- 
persion are not considered. Furthermore, Statler uses data only 
outside 15". However, Statler has obtained new, high S/N data 
(Statler & Smecker-Hane 1999), which he will present in a fu- 
ture paper. 

The studies above use low S/N, sparse kinematic data, or no 
kinematic data at all. In this paper, we report results from obser- 
vations of the full line-of-sight velocity distribution (LOSVD) 
at various position angles from high S/N ground-based data 
and data from HST. The data are fit to axisymmetric, three- 
integral orbit-based models. Richstone et al. (1999) describe 
the modeling technique in a companion paper. We exploit the 
full LOSVDs directly in the modeling. Alternatively, we could 
use the moments of the LOSVDs; however, the full shape of 
the LOSVD contains important information on the underlying 
dynamics, and it is desirable to incorporate the full LOSVD 
when the spectra have sufficient signal-to-noise to determine 
this reliably. In fact, our result for NGC 3379 critically de- 
pends on the shape of the central LOSVD. Cretton et al. (1999) 
and van der Marel et al. (1998) present a very similar model- 
ing technique, where the only difference from ours is that we 
fit the full LOSVD as opposed to their moment fitting proce- 
dure. There are other programs that also use the higher mo- 
ments of the LOSVD as constraints (Dejonghe, 1987; Gerhard 
1993; Rix et al. 1997). Obtaining the full LOSVD at many 
position angles is the maximum obtainable kinematic informa- 
tion and has now become routine {e.g. Bender et al. 1994, Car- 
olio et al. 1995, Gebhardt & Richstone 1999). 

In §2, we describe the photometric and kinematic data for 
NGC 3379. In §3, we describe the models and present results 
in §4. Discussion is given in §5. 



2. DATA 

2.1. Photometry: HST and Ground-Based Imaging 

We obtained Hubble Space Telescope (HST) observations 
of NGC 3379, consisting of four 500s F555W (V-band) and 
four 400s F814W (/-band) WFPC2 images, on 1994 Novem- 
ber 19. Centering the galaxy in the high-resolution PCI camera 
(0.0455" pixel" 1 ) provided a total signal in the central pixel of 
~2x 10 4 e~ in both V and /. We deconvolved the summed im- 
ages with 40 Lucy-Richardson iterations. The surface bright- 
ness, ellipticity, and position angle of the major axis were ob- 
tained using the same procedure as in Lauer et al. (1995) and 
are presented in Table 1 . 
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FIG. 1 . — Visual surface brightness and luminosity density for NGC 3379. 
The solid circles are the data from HST and the solid triangles are the data from 
Peletier et al. (1990), which have been adjusted to match with the HST data 
assuming constant color. The solid line in the left panel represents a smooth- 
ing spline (Gebhardt et al. 1996) of the surface brightness data. The inferred 
luminosity density is given for three different inclinations; from bottom to top 
the lines represent the densities at 90° (intrinsic El), 31° (E4), and 27° (E6). 

The HST profile reaches only about 12" from the center and 
must be extended to larger radii using ground-based photom- 
etry. Peletier et al. (1990) have obtained ground-based UBR 
photometry for NGC 3379 over the radius range 4"-170", and 
Capaccioli et al. (1990) have B photometry out to 500"; we 
must therefore estimate the transformations to V and / using the 
overlap region with the HST data. The data from Peletier et al. 
and Capaccioli et al. agree well, and since Peletier et al. provide 
several colors, we use their results. The colors in the overlap re- 
gion, based on the average differences there, are B-V = -0.58 
and R-I = -0.94. We have ignored possible variations in the 
colors as a function of radius. Using these colors, we com- 
bine the HST and ground-based data to provide the full profile. 
The left panel in Fig. 1 shows the surface brightness profile in 
V along the major axis from 0.02"-170". The three panels in 
Fig. 2 show the ellipticity, the position angle of the major axis 
(measured from north through east), and the color V-/ as func- 
tions of radius along the major axis. The random uncertainties 
can be estimated from the local scatter. 

The dynamical modeling described below assumes axisym- 
metry. We must therefore determine the best axisymmetric den- 
sity distribution consistent with the observed surface brightness 
distribution. Strictly speaking, there is none, because, as shown 
in Figure 2, there are abrupt changes in the ellipticity and po- 
sition angle at radii less than 1" and larger than 100". PA 
changes in particular are normally attributed to variation in the 
axis ratios with radius in a triaxial distribution, or to real three- 
dimensional twists in locally axisymmetric objects in which 
surfaces of constant density, while axisymmetric, are not co- 
axial. 
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FIG. 2. — Ellipticity, position angle and V — I color as a function of radius 
along the major axis for NGC 3379. The solid circles are the data from HST, 
and the solid triangles are data from Peletier et al. (1990). The ground-based 
V-I is estimated from B-R. 

However, for several reasons we believe that, at radii less 
than 100", an axisymmetric model will suffice for ourpurposes. 
First, in the middle range of radii, 1" < r < 100", the ellipticity 
and PA are fairly constant at 0.1 and 72 ± 5°. Second, in the 
central region there is a dust ring with radius 1" (discussed be- 
low and seen in Fig. 5). The absorption by the dust ring - 1% 
of the galaxy surface brightness - is enough to affect estimates 
of both the ellipticity and the position angle. Furthermore, the 
shallow central slope (dlog S/dlog R = -0.2) of the inner pro- 
file and the small ellipticity (<~ 0.1) render the position angle 
at small radii (R < 1") highly uncertain. Even if the PA gra- 
dient seen in Fig 2 is real and well determined, the amount of 
light, and therefore the gravitational impact, of the associated 
departure from axisymmetry is small. Finally, at large radii, 
departures from axisymmetry have little impact on the central 
gravitational field. 

The data thus suggest that NGC 3379 can be reasonably well 
represented by a spheroid of constant ellipticity. The deprojec- 
tion of the surface brightness is then unique for a given inclina- 
tion. The deprojection is based on a non-parametric estimate of 
the density using smoothing splines (Gebhardt et al. 1996). The 
right panel of Fig. 1 plots the luminosity density for three dif- 
ferent inclinations (90° is edge-on). For different inclinations, 
and hence different flattenings, the luminosity density scales as 
the ratio of the apparent to the true axis ratio. 

We did not examine models in which the galaxy is assumed 
to be axisymmetric but not spheroidal (i.e., with varying el- 
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lipticity), in part because in this case the deprojection is not 
unique. Romanowsky & Kochanek (1997) determine the uncer- 
tainties in the deprojection of an axisymmetric system for vari- 
ous inclinations, and the uncertainties can be quite large for low 
inclinations (close to face-on). We could test the effect of our 
spheroidal assumption by modeling the most extreme situation: 
a galaxy which is seen nearly face-on and containing an embed- 
ded disk (similar to the suggestion of Capaccioli et al. 1991). 
The most likely effect is to change to inclination and M/Ly con- 
straints, and relaxing the spheroidal assumption should have lit- 
tle effect on the measured black hole mass. This issue, however, 
will not be addressed here. 

2.2. HST FOS Spectroscopy 

NGC 3379 was observed on 1995, January 23 - 27 with the 
0.21" square aperture ("0 . 2 5-PAIR") of the FOS. The spec- 
trum was taken in two visits; the total integration time was 
215.5 minutes. 

The wavelength range, 4566-6815 A, includes the Mg I b 
lines at A ~ 5175 A and the Na I D lines at A ~ 5892 A. The 
spectrum was electronically quarter-stepped, giving 2064 1/4- 
diode pixels. The reciprocal dispersion was 1.09 A pixel" 1 . 
The instrumental velocity dispersion was measured to be 
cr instr =FWHM/2.35= 1.76 ± 0.03 pixels = 101 ±2 km s" 1 (in- 
ternal error). This width is intrinsic to the instrument and 
is not strongly affected by how the aperture is illuminated 
(Keyes et al. 1995). We therefore make no aperture illumina- 
tion corrections to the measured velocity dispersion. 

Flat fielding and correction for geomagnetically induced mo- 
tions (GIM) were done as in Kormendy et al. (1996). The 
flat-field image was taken with the same . 2 5-PAIR aper- 
ture used for NGC 3379. The galaxy exposure was divided 
into 3 and 4 subintegrations for the two visits. Because GIM 
shifts between subintegrations are small, unlike the large shifts 
between galaxy exposures and the flat field, we first averaged 
all the exposures and used cross-correlation to determine the 
shift between this average and the fiat field. The measured 
shift was 2.27 pixel. Each subintegration was multiplied by 
the fiat-field frame shifted by the above amount. One dead 
diode was corrected by dividing the affected pixel intensities 
by 0.8. Then GIM shifts between subintegrations were deter- 
mined; most shifts were < 0.5 pixels, with the largest being 0.8 
pixels. Each subintegration was then shifted to agree with the 
subintegration taken closest in time to the comparison spectrum 
exposure. Finally, the spectra were added, weighted by the ex- 
posure times. 

2.2.1. Extracting the LOSVD 

Obtaining the internal kinematic information requires a de- 
convolution of the observed galaxy spectrum with a represen- 
tative set of template stellar spectra. We deconvolve the spec- 
trum using a maximum-penalized likelihood (MPL) estimate 
that obtains a non-parametric line-of-sight velocity distribu- 
tion (LOSVD). The fitting technique proceeds as follows: we 
choose an initial velocity profile in bins. We convolve this pro- 
file with a weighted-averaged template (discussed below) and 
calculate the residuals to the galaxy spectrum. The program 
varies the velocity profile parameters — bin heights — and the 
template weights to provide the best match to the galaxy spec- 
trum. The MPL technique is similar to that used in Saha & 
Williams (1994), and Merritt (1997). However, in constrast 
to these authors we simultaneously fit the velocity profile and 



template weights. We impose smoothness on the LOSVD by 
the addition to the \ 2 of a penalty function, which is the in- 
tegrated squared second derivative (details are in Merritt 1997 
and Gebhardt & Richstone 1999). We estimate the best smooth- 
ing from bootstrap simulations (described below). Two tem- 
plates, DSA107-442 (a field K0III star) and N188-I69 (a K2III 
star in NGC 188), are used, both observed with the FOS. Star 
N 188-169 is the preferred template, as the fitting routine gives 
essentially all of the weight to this star. However, the LOSVD 
determined using only the other star is similar, suggesting that 
the measurement of the LOSVD is robust. The results from this 
maximum likelihood technique were compared to results of a 
Fourier correlation technique (Bender 1990), and no significant 
differences were found. 

The LOSVDs obtained in the two visits by HST agree to 
within their uncertainties, and also the pointings from the two 
visits agree to within 0.01". Therefore, we sum the two spectra 
weighted by exposure time to derive a combined LOSVD. A 
parameterization of the LOSVD by Gauss-Hermite moments 
(van der Marel & Franx 1993, Gerhard 1993) gives for the 
first four moments v = ± 14 km s" 1 , a = 289 ±13 km s" 1 , 
h 3 = 0.1 1 ± 0.04, and h A = 0.04 ± 0.04. The velocity zero-point 
is arbitrary and will be discussed in §4. All error bars corre- 
spond to the 68% confidence band. 

In the lower part of Figure 3 we plot the template spec- 
trum; the noisy upper line is the averaged FOS spectrum for 
NGC 3379 (shifted by its redshift), and the smooth upper line is 
the template spectrum convolved with the derived LOSVD. The 
continuum has been divided out of both template and galaxy 
spectra. The spectrum shown here is a region around the Mg I b 
triplet (5175 A) and not the full spectrum obtained by the FOS. 
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FIG. 3. — Spectra for the template star (lower line), the central 0.21" of 
NGC 3379 (noisy upper line), and the template star convolved with the de- 
rived central LOSVD (smooth upper line). 

The 68% confidence band for the LOSVD is determined 
using a bootstrap approach. We convolve the template star 
with the measured LOSVD to provide the initial galaxy spec- 
trum (the upper smooth line in Fig. 3); from that initial galaxy 
spectrum, we then generate 100 realizations and determine the 
LOSVD each time. Each realization contains randomly chosen 
flux values at each wavelength position drawn from a Gaussian 
distribution, with the mean given by the initial galaxy spectrum 
and the standard deviation given by the root-mean-square of 
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the initial fit. The average root-mean-square for the continuum- 
divided spectrum is 0.033. The 100 realizations of the LOS VDs 
provide a distribution of values for every LOSVD bin, from 
which we estimate the confidence bands. 
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FIG. 4.— LOSVD for NGC 3379. The solid and dotted lines are the best 
estimate and the 68% confidence band from the central HST point. The dashed 
line is a symmetrized LOSVD. 

2.2.2. Asymmetry in the LOSVD 

Figure 4 presents the LOSVD with the 68% confidence band 
(solid and dotted lines). The LOSVD is significantly skewed 
with a tail to positive velocities - a feature also detected using 
the Fourier correlation technique. The skewness is crucial for 
the following analysis and must be considered carefully. In an 
axisymmetric and even a triaxial system, the central LOSVD is 
symmetric. 

Various attempts were made to get rid of the skewness by 
reducing the data differently; we used both templates individ- 
ually, we used flattened and un-flattened spectra, we used the 
Mg b lines and the Na doublet separately and combined, dif- 
ferent shifts were used to combine the spectra, and the spectra 
from both peak-ups were reduced separately. In every case, the 
skewness remained. Since the skewness appears to be real, we 
must conclude that the FOS was not centered, that part of the 
galaxy light is being obscured, that the galaxy is not symmetric 
near the center (possibly because of an offset nucleus), or any 
combination of these affects. If the nucleus is offset, we can- 
not successfully model NGC 3379 as an axisymmetric system. 
However, there is no evidence to suggest that NGC 3379 has an 
offset nucleus since the isophotes, except for the dust ring, are 
concentric. 

Figure 5 is an image of the central 5" of the galaxy with 
the spheroidal component subtracted. We include in Fig. 5 the 
exact position of the FOS aperture. To determine the exact po- 
sition of the aperture, we have two independent checks: from 
the values used to peak-up on the galaxy center, and from the 
acquisition image taken immediately before the galaxy obser- 
vations. For the peak-up procedure, HST steps by 0.065" and 
points the FOS at the highest pixel value. Examination of the 
peak-up numbers reveals that the likely miscentering is 0.033", 
since the center of the peak-up numbers lies roughly between 
two bins in the peak-up array. This miscentering is also con- 
firmed from the acquisition image which places the FOS aper- 
ture the same distance from the galaxy center. 
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FIG. 5. — Image of NGC 3379 with the spheroidal light distribution sub- 
tracted. The dust ring is clearly visible. The major axis of the galaxy is shown 
as the solid line. The box is the size of the FOS aperture (0.21" square). The 
uncertainty in the position of the box is approximately one pixel (0.045"). The 
scale ranges from -200 to +100 electrons per pixel. 

In order to explain the LOSVD skewness, the direction of the 
aperture offset must correspond appropriately. HST placed the 
aperture on the eastern side of the galaxy, along the major axis. 
Thus we acquire more light from the eastern side relative to the 
western side. The NE side of the galaxy has negative rotation 
velocities (Davies & Birkinshaw 1988). If one side contributes 
more light and there is significant rotation, the mode and mean 
of the distribution shifts towards that side — in this case the NE 
side with negative velocity — and the other side still contributes 
an extensive tail toward positive velocities, thus generating the 
observed skewness in the LOSVD. 

A further source for the LOSVD skewness may be that part 
of the galaxy light is obscured. A dust ring is clearly visible, 
inclined by about 45° from the major axis (the solid line in 
Fig. 5). From the relative depths of the obscuration on oppo- 
site sides, it appears that the SW side of the ring is on the near 
side of the galaxy, while the NE side is on the far side. This 
configuration means that the crucial signal from stars near the 
nucleus is preferentially obscured on the SW side, thus leading 
to the same effect as the miscentering; the western side con- 
tributes less light than the eastern side causing positive velocity 
skewness. 

The miscentering of the aperture results from two effects, 
errors in the peak-up procedure and errors due to the dust ob- 
scuration. To calculate the miscentering from dust, we find the 
best center by using the galaxy light beyond the dust ring and 
compare this with the FOS center. HST points FOS to the high- 
est position in a 0.21" aperture on a 4 x 4 raster with spacing of 
0.065". The dust does not appear to have a significant effect on 
the location of the brightest pixel; it causes a miscentering of 
only 0.01". However, errors in the peak-up procedure are more 
significant, and as mentioned above, the miscentering can be 
0.033". We conclude that the maximum miscentering is 0.04" 
(peak-up errors + dust miscentering). We then place the FOS 
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aperture in the WFPC2 image at a position 0.04" East from the 
galaxy center and determine the relative amount of light in the 
aperture from opposite sides of the galaxy. At this position, the 
light from the NE side of the galaxy is twice as strong as from 
the SW side. Thus, in the central bin of our model, we include 
all of the light from the NE side and only half from the S W side. 

An alternative is to use the best symmetrized LOSVD in the 
modeling and assume equal contributions from both sides of 
the galaxy. The dashed line in Fig. 4 is the maximum penal- 
ized likelihood estimate of the best symmetric LOSVD. Results 
from both assumptions are discussed below. 

2.3. Ground-Based Spectroscopy 

Ground-based spectroscopic data are available from Geb- 
hardt & Richstone (1999) (along four position angles), 
Statler & Smecker-Hane (1999) (along four position angles), 
Bender et al. (1994) (major axis only), Gonzalez (1994) (major 
and minor axes), Franx et al. (1989) (major and minor axes), 
and Davies & Birkinshaw (1988) (major, minor, and 30° axes). 
The data from Franx et al. and Davies & Birkinshaw have lower 
S/N and will not be considered here. The data from the first four 
sources all have high signal-to-noise (S/N) and in general agree 
with one another (see Gebhardt & Richstone for a detailed com- 
parison), with the largest disagreement coming from the large 
radii data. For the data of Gebhardt & Richstone, the position 
angles observed and the radial extraction of the galaxy spectra 
were designed to correspond exactly with the binning used in 
the modeling (described in the next section), and we use their 
data. 

3. MODELS 

We use 3-integral axisymmetric models (Richstone et al. 
1998). The models are orbit-based (Schwarzschild 1979) and 
constructed using maximum entropy (Richstone & Tremaine 
1988). We run a representative set of orbits in a specified po- 
tential, and determine the non-negative weights of the orbits 
to best fit the available data. Maximizing the entropy helps to 
provide a smooth phase-space distribution function. Rix et al. 

(1997) present a similar code which has been applied to the 
spherical case; Cretton et al. (1999) and van der Marel et al. 

(1998) developed a fully-general axisymmetric code that is also 
orbit-based. The difference between our code and theirs is due 
to how we impose smoothness and how we fit to observational 
data. Smoothness in our code comes from maximizing entropy, 
whereas the other authors enforce smoothness of the distribu- 
tion function directly by minimizing its variation. The most im- 
portant difference is that we fit the LOSVD directly, and both 
Rix et al. and Cretton et al. fit moments of the velocity profile. 
Using the LOSVD precludes the need to approximate the veloc- 
ity profile with a parametric function both in the observations 
and in the orbital libraries. 

3.1. Generating Orbit Libraries 

The surface brightness yields an estimate of the luminosity 
density distribution as explained above (§2.1), assuming that 
the luminosity density is axisymmetric on similar spheroids and 
that the galaxy has a given inclination. To obtain the mass den- 
sity distribution, we have made one of two assumptions about 
the mass-to-light ratio (M/Ly): (i) constant M/Ly, (ii)M/Ly that 
varies according to the V -I color shown in Fig. 2. NGC 3379 
is redder in the center by A(V -I) ~ 0. 1 . The analysis of Gon- 
zalez (1994) indicates that this color change is due mainly to 



metallicity (i.e., no age effects). In this case, models of Worthey 
(1994) predict that M/Ly may be as much as 25% higher in the 
center. The V-I variation in Fig. 2 is approximately log-linear, 
so we vary M/Ly by 25% from the outer to inner radii following 
the same functional form. Finally, we add a central point mass 
and compute the potential. 

Using the derived potential, we follow orbits which sample 
the available phase space in energy (E), angular momentum 
(L z ), and the third integral (7 3 ) (Richstone et al. 1998). Each 
orbit crosses the equatorial plane 40 times. Binning the orbital 
distributions both in real space and velocity, and in projected 
distance and line-of-sight velocity, provides an estimate of the 
contribution to the LOSVD from each orbit at each spatial bin. 
Both the real and projected spatial binning are approximately 
logarithmic in radius and linear in cos 9, where 9 is measured 
from the pole. This binning scheme provides approximately 
equal-mass bins. Since we are using axisymmetric modeling, 
we follow orbits only in r and 9. We therefore have to choose 
the azimuthal angle for the projection, and, since all angles are 
equally likely, we choose that angle at random. For every time 
step, 100 azimuthal angles are chosen. 

We measure an orbit's contribution on a grid of 80 radial 
bins, 20 angular bins, and 13 velocity bins — compressed to 
20 x 5 x 13 when comparing to the data. The size of the central 
bin equals the aperture used in the FOS (0.21"). The outermost 
bin ranges from 149" to 200". The centers of the five angular 
bins up from the major axis are at 6°, 18°, 30°, 45°, and 72°. 
The velocity bins are about 100 km s" 1 wide, approximately 
equal to the spectral resolution of the ground-based observa- 
tions. 

We follow <~ 3200 orbits per model depending on inclination 
and black hole mass. This number is only for one sign of the 
angular momentum; consequently, it must be doubled before 
comparison to the data. We flip the LOSVD at every position 
in the galaxy to provide the appropriate reversed orbit. Thus, 
our final model contains ~ 6400 orbits. Varying the relative 
weights of these orbits allows us to determine the best match to 
the data. 
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3.2. Modeling the Data 

It is desirable to constrain every model bin with an observed 
LOSVD. For NGC 3379, we have ground-based observations 
of the LOSVD along four of the five position angles used in the 
modeling. In fact, the observations were designed to mimic the 
model angular bins listed above. Furthermore, the spectra were 
extracted from the data using the same radial binning scheme as 
the modeling. Therefore, most model bins correspond exactly 
to an observation; i.e., no interpolation of the data or the model 
is necessary for comparison. This correspondence breaks down 
for the ground-based data in the innermost bins because of see- 
ing. 

The ground-based data extend out to 64", about two effective 
radii. The total number of bins for which we have an observa- 
tion is 54. The model LOSVDs are calculated at 13 velocity 
positions with a spacing of 100 km s -1 , providing at most 13 
independent measurements of the LOSVD per position. We 
therefore have a maximum of 702 (54 positions x 13 LOSVD 
bins) independent observables. However, the actual number 
is smaller than this, mainly due to the smoothing used in the 
measurement of the LOSVDs (discussed in §5 and Gebhardt & 
Richstone 1998), which correlates the 13 velocity bins. Thus, 
we will overestimate the goodness-of-fit since we use the 68% 
confidence bands of the LOSVDs to calculate it. 

The model incorporates seeing by convolving the light dis- 
tribution for every orbit with the appropriate PSF before com- 
paring to the data. For the ground-based data, the estimated 
PSF, approximated as Gaussian, has a FWHM of 1.5", includ- 
ing both atmospheric seeing and slit size. The HST data point 
was taken in a square aperture of 0.21". The PSF for HST has 
a FWHM of around 0.07", and hence smearing due to the PSF 
is negligible compared to the aperture used. Therefore, for the 
HST data, only the aperture size was considered with no PSF 
convolution. We have checked in a number of models that this 
approximation has no effect. 

For each model, we match both the light in each of the 20 x 5 
real-space bins, and the LOSVDs in the 54 locations where we 
have data. The technique minimizes the \ 2 between the model 
and data LOSVDs. The light from the model matches the de- 
convolved spatial brightness profile to a tolerance of better than 
1 % in each bin, which also ensures that the surface brightness 
matches better than 1 % in any projected light bin. The match 
to the total light in each bin has to be done using the real spatial 
profile and not the surface brightness profile. If the matching 
were done in projection, it would admit different density distri- 
butions than the one for which the orbits were computed since 
the deprojection of a flattened, inclined, axisymmetric system 
is not unique. 

4. RESULTS 

4. 1 . Estimating the Goodness-of-Fit 

For each model, the 702 observables (13 LOSVD bins at 
the 54 positions) determine the total % 2 , given by X^LitCy* - 
yi)l Ci) 2 where the y's are the LOSVD bin heights for the model 
and the data and a is the uncertainty from the 68% confidence 
band. As an example, Figure 6 plots the LOSVDs for both the 
data and the best model at four positions: the center, a major- 
axis position near to the center, a position close to one effective 
radius on the major axis, and a minor-axis position near to the 
center. 

Judging the goodness-of-fit from the \ 2 values of the various 



models is not straightforward because the number of the de- 
grees of freedom is not easy to determine; the derivation of the 
LOSVDs includes a smoothing parameter which correlates the 
values and uncertainties of the velocity bin heights (Gebhardt & 
Richstone 1999). With higher S/N data than present, we could 
reduce the smoothing parameter to lessen this problem. We do 
not attempt to estimate directly the actual degrees of freedom 
and therefore do not have an overall goodness-of-fit measure, 
but instead we calculate the change in \ 2 as a function of the 
three variables - black-hole mass, inclination, and M/Ly. The 
lowest x 2 value thus provides the best model, with the uncer- 
tainties given by the classical estimators for the A% 2 's. 
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FIG. 6. — Projected LOSVDs at four positions: the central HST/FOS posi- 
tion (top-left), along the major axis near the center (top-right), along the major 
axis near an effective radius (bottom-left), and along the minor axis near the 
center (bottom-right). The data are the open circles with their corresponding 
eiTor bars. The solid points are the model values. The area is normalized to 
the total light in that bin. This model is our preferred model, with 1 X 10 s Mq 
black hole and inclination 90° . Note that the central HST LOSVD has been 
nipped about zero relative to Fig. 4 in order to measure the model with respect 
to positive rotational velocity. 



The total x 2 's for most of the models have a value near 200, 
giving a reduced \ 2 of 0.3, assuming one observable equals one 
degree of freedom. A reduced \ 2 near unity would require ei- 
ther that every three LOSVD bins are correlated by the smooth- 
ing, therefore reducing the number of degrees of freedom by 
three, or that the uncertainties of the bin heights are overesti- 
mated by a factor of y/3. In the former case, using Ax 2 as a 
discriminator for the error of the best-fitting models will not 
be affected since the number of degrees of freedom has no in- 
fluence on the calculation of the \ 2 differences; A\ 2 = 1 will 
still correspond to 1-er. However, in the latter case, both the % 2 
and Ax 2 values would increase by a factor of three, creating 
a smaller confidence band on the variables; A\ 2 =1/3 would 
then correspond to 1-er. It is likely, however, that the dominant 
reason for the low reduced \ 2 lS the underestimation of the de- 
grees of freedom since we know a priori that the smoothing 
must have an effect. To provide a quick check on smoothing 
effects, we refit the LOSVDs for a subset of positions using a 
much smaller smoothing length. As expected, the LOSVD be- 
comes noisier. Measuring the \ 2 between this new LOSVD and 
the model that we fit using the the older LOSVD, we find that 
the increase in \ 2 ranges from 3-5. Thus, we conclude that the 
dominant effect of the smoothing is to correlate the variables 
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and reduce the number of degrees of freedom. In retrospect, it 
would have been better not to have smoothed the LOSVDs and 
preserve their statistical independence. Note however that our 
decision to accept Ax 2 = 1 as the 1-cr error is conservative and, 
if anything, overestimates the errors. 



4.2. Model Fits 

We examined models with inclinations from 29° to 90° (in- 
trinsic E5 to El galaxy), black holes from zero to 10 9 M Q , and 
M/Ly from 4.0 to 8.0 in solar units. The x 2 is a function of these 
three variables. For each inclination, Fig. 7 plots x 2 as a func- 
tion of black-hole mass (zero black-hole mass is shown as 10 6 ), 
using that M/Ly which provides the smallest x 2 at mat black- 
hole mass. All inclinations show a minimum x 2 a t black-hole 
masses in the range 6 x 10 7 -2 x 10 8 M Q . The fits clearly prefer 
the edge-on models (90°). Fig. 8 presents contours of A% 2 as a 
function of black-hole mass and M/Ly for each inclination. The 
points represent locations of the modeled values. The contours 
use a two dimensional smoothing spline (Wahba 1980) to esti- 
mate the x 2 values at positions where no models occur. As in 
Wahba, Generalized Cross-Validation determines the smooth- 
ing value; however, the modeled values are relatively smooth 
and little smoothing is necessary. 
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FIG. 7. — x 2 as function of black hole mass for various inclinations. The la- 
bels give the inclination and the intrinsic galaxy shape. We have added 10 6 Mq 
to the models with no black hole to put them on the log scale. 



Since we have such a large parameter space — minimizing 
several thousand variables — we must ensure that the model- 
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FIG. 8. — Contours of A\ 2 as function of black-hole mass and M/Ly for four inclinations. The solid points represent the models that were run. The contour 
values are equal to A% 2 = 1,2. 7, 4.0, and 6.6; one-dimensional projection of these contours provide the 68%, 90%, 95% and 99% confidence bands respectively for 
the two variables, black-hole mass and M/Ly . 
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ing program is finding a true minimum. Furthermore, we re- 
quire that the numeric noise caused by the minimizer be smaller 
than the quoted uncertainties of the output parameters (e.g., the 
black hole mass). The two conditions that govern these checks 
are the sampling densities that we use in parameter space, and 
the tolerance used to determine convergence in the minimizer. 
Fig. 8 demonstrates the sampling densities for the black-hole 
mass and the M/Ly, where each point represents one model. 
Fig. 7 shows the range and number of modeled inclinations. 
The smoothness of the \ 2 contours in both figures demonstrate 
that we have found the global minimum. As a further check, 
we have run models with both larger and smaller M/Ly than 
those presented, however the fits are significantly worse and 
we do not include them in the plots. Thus, we conclude that 
the sampling density provides adequate coverage around the 
global minimum. Our second concern is whether the tolerance 
used for the stopping criteria in the minimizer creates signifi- 
cant noise for the parameter estimation. We must use a stop- 
ping tolerance in the program or the time it takes for the x 2 
to asymptote becomes impractically long. Both Figures 7 and 
8 provide an estimate of the numeric noise. From the sooth- 
ness of plots, the change in x 2 between two adjacent points is 
much smaller than the global change in % 2 . Also, in a handful 
of models, we iterate till there is no change in x 2 relative to 
machine precision, and find that this asymptotic value is within 
Ax 2 = 0.02 of the stopping value. Therefore, the noise in x 2 
generated from the minimization routine is not significant. The 
smoothness is the solution is also seen in Fig. 9, where we plot 
best-fit black hole mass as a function of inclination. Although, 
the sampling of points is not dense, the fact that there are no 
abrupt changes signifies we are not subject to numeric noise in 
the solution. 
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FIG. 9. — The best-fit black-hole mass as a function of inclination. 





r (arcsec) 

FIG. 10. — Internal dynamics of NGC 3379 for the preferred model, 
1 X 10 s Mq black hole, inclination 90°. The top plot shows the three internal 
dispersions along the major axis; the bottom plot shows the mean velocity in 
the equatorial plane. 

As described in §3.1, we have also run models with M/Ly in- 
creasing 25% log-linearly with radius towards the galaxy cen- 
ter. In this case, the required M/Ly is systematically lower than 
when using the constant M/Ly assumption; however, the results 
for the required black hole and inclination are unchanged. 

The best model has an inclination of 90°, a black hole mass 
of 1 x 10 8 M Q , and M/Ly= 4.85. The allowed mass of black 
holes range from 6 x 10 7 -2 x 1O 8 M . The internal dynam- 
ics of the best-fitting model are shown in Fig. 10. The top 
plot shows the internal velocity dispersions of the three com- 
ponents along the major axis. The model is tangentially bi- 
ased (in the (f> direction) in the smallest bin, radially biased in 
the mid-range radii, and tangentially biased again at the largest 
radii. Independent of inclination, Fig. 10 well characterizes the 
orbital distribution for each of the best-fitting models. Simi- 
larly, the mean velocity in the equatorial plane (bottom plot of 
Fig. 10) has a general shape for the best models; it shows no 
systematic trends from 0.4-30" (20-1500 pc), increasing sig- 
nificantly in the center due to the presence of the massive black 
hole. Beyond about 60" the model is unconstrained since there 
is no kinematic information, but we plot the model results there 
since they represent the maximum entropy configuration. For 
each inclination the trend with black-hole mass is as follows: 
increasing the black-hole mass requires an increased amount of 
tangential anisotropy in the central bin, a decrease of the radial 
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FIG. 11. — Gauss-Hermite moments for data — the points with 68% error bars — and models. The solid line is the best-fit edge-on black-hole model, the dased blue 
line is the edge-on no-black-hole model, and the dotted orange line is the best-fit 29° model. Each column represents a different position angle which is measured 
up from the major axis. 



bias in the mid-range radii, and the orbits at the largest radii 
to become more isotropic. The coupling between the black- 
hole mass and the amount of radial anisotropy in the mid-radii 
stems from the need to match the same central LOSVD; a large 
central dispersion can be modeled with a large black-hole mass 
or radial orbits. Therefore, as the black-hole mass increases, 
the need for radial orbits decreases. In the very central bin, 
a large black-hole mass induces large rotational velocities and 
also large dispersion in the <fi direction. At the largest radii, the 
situation is more complex since we have not included a dark 
halo; if the projected dispersion is larger than the dispersion ob- 
tained assuming isotropy and the stellar mass profile then one 
can either invoke a dark halo or include tangential anisotropy. 
In any event, it is not clear whether meaningful analysis at large 
radii can be obtained from models which do not include a dark 
halo. 

To demonstrate why the data seem to prefer an edge-on 
model that contains a black hole, we must look at the model 
fits to the data. Figure 11 plots Gauss-Hermite polynomial 
estimates for the data and model LOSVDs for all four posi- 
tion angles. The difference between the edge-on model and the 
nearly face-on model is almost entirely at large radii; the edge- 
on model matches the large-radial data better, specifically for 



the velocity and velocity dispersion. In contrast, the difference 
between the no-black-hole and black-hole model is so subtle 
that one can barely discriminate those two models in Fig. 11. 
Only in the central bin does the the black-hole model slightly 
better estimate the skewness in the profile. However, Fig. 1 1 is 
not the optimal way to compare models since the fit uses the full 
LOSVD and not moments. Instead of plotting all 54 LOSVDs 
for various models, Fig. 12 plots the difference in \ 2 at each 
spatial position. This difference represents the \ 2 measured 
for the LOSVD of a particular model minus that of the best-fit 
model (edge-on, black-hole model). The two comparison mod- 
els are the same as in Fig. 11, the best no-black-hole model (an 
edge-on model) and a nearly face-on black-hole model (at 29°). 
As expected, the discrimination between black-hole/no-black- 
hole comes mainly from the HST/FOS data point. A no-black- 
hole model fails to match the shape of the central LOSVD; how- 
ever, there is a limit to size of the central black hole since very 
high black-hole masses models over shoot the dispersions in 
the central few ground-based measurements. Without the HST 
data point, we would not have been able to determine whether 
a black hole exists in NGC 3379. The discrimination for the in- 
clination comes from the data at larger radii; here, the edge-on 
model is much better able to match the projected kinematics. 
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As we discuss later, this inclination constraint is uncertain due 
to our assumption of constant M/L v . 
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FIG. 12. — The change in \ 2 relative to the best-fitted model. Each panel 
represents a different observed position angle. The solid line represents the \ 2 
difference for the edge-on no-black-hole model, and the dashed line represents 
the 29° model with a 1O 8 M black hole. 

There is an uncertainty in the absolute wavelength calibra- 
tion of the HST spectrum. Since we do not assume that the 
central HST LOSVD is symmetric, we must determine the ve- 
locity zeropoint of this profile relative to the rest frame of the 
galaxy (920± 10km s" 1 , Borne et al. 1994). Unfortunately, our 
original observing strategy did not include an accurate wave- 
length calibration, since we were concerned mainly with the 
central dispersion, and thus did not obtain the necessary flank- 
ing arc-lamp spectra. We therefore tried a range of velocity 
zeropoints, models with velocity offsets from 0-50 km s" 1 . 
As the offset increases, the lower limit on the black-hole mass 
grows, whereas the upper limit remains the same. We are not 
able to fit acceptable models with offsets greater than 15 km s -1 . 
Therefore, our allowed range of black hole masses is largest us- 
ing the km s" 1 velocity offset. All of the figures contain the 
results using this value. 

5. DISCUSSION 

The models in this paper depend primarily on the input black- 
hole mass, inclination, and M/L assumption. A difficulty is that 
different combinations of these variables can have similar ap- 
pearances in the projected kinematics. For example, the incli- 
nation and the M/L variation both affect the fit at large radii; one 
can easily see how differently inclined models of the same flat- 
tened galaxy will have markedly different projected dispersion 
profiles, and also the dispersion profile is obviously influenced 
by the inclusion and shape of a dark halo. Thus, by varying the 
shape and amount of dark halo one may be able to mimic the 
effect of changing the inclination. Since we have not explored 
a full range of halo models, we will not make strong claims as 
to our ability to constrain the inclination of the galaxy. How- 
ever, the black hole mass is mainly driven by the data at small 
radii, and primarily by the HST/FOS data point. Fig. 9 demon- 
strates this independence from the inclination, as our best-fit 
black hole mass does not depend strongly on the inclination. 
An even more robust measurement is the orbital structure. As 
either the inclination or the black-hole mass varies, we find a 



similar velocity ellipsoid; at the smallest bin the velocity ellip- 
soid is mainly tangential, becoming radial in mid-range radii, 
and tangential again at the largest radii. Fig. 13 plots the radial 
motion relative to the tangential motion for four models at dif- 
ferent inclinations. Determination of the inclination, black-hole 
mass, and orbital structure is discussed in detail below. 
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FIG. 13. — Radial relative to tangential dispersion as a function of radius 
along the major axis for the best-fit models at the specified inclination. The 



tangential dispersion is defined as <r 2 



= («2+«i)/2. 



The inclination of NGC 3379 has been the subject of some 
debate. Capaccioli et al. (1991) conclude that NGC 3379 is 
similar to NGC 3115 but nearly face-on, with some degree of 
triaxiality. Capaccioli base this conclusion only on the sur- 
face brightness profile and use no kinematic data. In contrast, 
Statler (1994) concludes that nearly face-on, triaxial models are 
inconsistent with the data at the 98% confidence limit, strongly 
contradicting Capaccioli's result. He prefers less-inclined mod- 
els. Statler has obtained more recent data (Statler & Smecker- 
Hane, 1999) and a fuller analysis is in preparation. 

As seen in Sect. 4, our best models prefer an edge-on con- 
figuration. Preference for edge-on models arises from two ob- 
servations; first, the dispersion profiles along different axes are 
slightly different, and, in an axisymmetric case, this configura- 
tion is not possible for a face-on system. The difference mainly 
arises from the inability to match the large radial data, however 
the face-on model also cannot match the shoulder in the disper- 
sion profile at 10-20" (previously noted by Statler & Smecker- 
Hane, 1999); the shoulder is more pronounced in the central 
position angles, and the edge-one models match this configura- 
tion well (see Fig. 11). Second, major-axis rotation exists which 
is not allowed for a face-on configuration of an oblate axisym- 
metric system. Thus these two restrictions are mainly a result 
of our axisymmetry assumption. 

More general models can be constructed than those presented 
here. We could in principle extend the modeling technique to 
include triaxial distributions but then the parameter space be- 
comes too large for us to search (each axisymmetric model 
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takes approximately 20 hours on an Ultrasparc). Second, our 
assumption that the galaxy is spheroidal could be relaxed. Ro- 
manowsky & Kochanek (1997) present a technique which pro- 
vides a range of different density distributions which similarly 
project. Third, we have assumed that M/Ly is independent of 
radius (i.e., there is no dark halo). Ciardullo et al. (1993), 
using radial velocity measurements of planetary nebulae, find 
no evidence for a dark halo in NGC 3379 out to 3.5 effective 
radii (120"); however their uncertainties are large and could, 
in fact, allow for some dark halo (Tremblay et al. 1995). In 
all of our models, the last velocity dispersion measurement (at 
64") along each position angle is higher than the model ve- 
locity dispersion, possibly suggesting the need for an increase 
in M/Ly. At one effective radius (34"), there is no need for 
an increase in M/Ly, which is consistent with previous studies 
of other galaxies (Kormendy & Westphal 1989, Kormendy & 
Richstone 1992, Kormendy et al. 1997). Thus, assumptions 
of axisymmetry, spheroidal distribution, and constant M/Ly all 
lead to uncertainties in the best-fit inclination. Detailed mod- 
eling including dark halos has been carried out by Rix et al. 
(1998) for NGC 2434 and Saglia et al. (1999) for NGC 1399, 
but only as spherical systems. Only by allowing for axisymmet- 
ric or triaxial shapes with a dark halo can we set more realistic 
limits on inclination. 

Even though we are not able to place limits on the inclina- 
tion, we can measure black-hole limits. The presence of a dark 
halo will have no effect on the need for a central black hole, or 
on its mass. The need for the black hole is directly a result of 
the shape of the central LOSVD observed with HST The large 
dispersion and the skewness drive the models to require a black 
hole. We have also run models assuming a symmetrized cen- 
tral LOSVD (the dashed line shown in Fig. 4). The allowable 
range of black hole masses remains the same for the inclined 
models in this case. However, for the edge-on case the model 
with no black hole fits as well. The symmetrized LOSVD has a 
dispersion of 275 km s" 1 , compared with 289 km s" 1 for the ac- 
tual LOSVD, and has less tail weight (as can be seen in Fig. 4). 
These differences permit a no-black-hole model. The central 
bin most influences the model; most of the orbits (> 70% for in- 
clined models and > 95% for edge-on models) contribute light 
to this bin, and, furthermore, many have their pericenter there. 
Therefore, the characteristics of the central LOSVD, in partic- 
ular the high-velocity components, have a strong effect on or- 
bits throughout the galaxy. In the edge-on case with the sym- 
metrized LOSVD, the orbits redistribute themselves in such a 
way as to create a slightly larger central dispersion (275 km s" 1 ) 
without the need for a black hole, and at the same time main- 
tain a fairly flat projected dispersion profile; i.e., at small radii, 
the orbits are mainly radial, whereas at large radii they be- 
come primarily tangential. However, the actual LOSVD war- 
rants matching the skewness and the high-velocity tails. With 
no black hole present, the model must use stars at larger radii to 
create the high-velocity components, thereby limiting its ability 
to redistribute the orbits at large radii and match the kinematics 
there. With a black hole, the black hole provides much of the 
high-velocity tail of the LOSVD for orbits locally and gives the 
model more freedom to better match the larger-radii data. The 
symmetrized and skewed LOSVDs are inconsistent with each 
other with greater than 95% probability, both using an estimate 
from the binned values and using the Gauss-Hermite moments. 
Given the arguments in §2.2.1, we must include the skewness of 
the central LOSVD and thus do not consider the symmetrized 
LOSVD models further. 



Relaxing the other two assumptions — axisymmetry and 
spheroidal distribution-is unlikely to change the black-hole 
limits. The models with no black holes all are extremely poor 
fits to the data; therefore, even if we allow for the variation in 
deprojected densities given by Romanowsky & Kochanek, the 
results are unlikely to change. Of course, since the deprojection 
is unique for the edge-on models, our results are general in the 
edge-on case. Triaxial models will allow even greater freedom 
to match the data, however it is not clear whether triaxiality is 
a viable configuration in the centers of galaxies (Merritt& Frid- 
man 1996, Merritt 1999), and we do not consider these models. 

The position angle of the dust ring in Fig. 5 is 45° away from 
the major axis of the galaxy. A dust ring in an axisymmetric 
system with an offset position angle is only neutrally stable but 
can be stabilized by a small amount of triaxiality. If the galaxy 
is nearly face-on, then triaxiality is not necessary to stabilize 
the dust ring; in the face-on configuration, since the dust ring 
is close to edge-on (about 75°), it must be in an orbit close to 
polar, and thus can have arbitrary azimuthal angle. Velocity in- 
formation along the dust ring could provide an important check 
on the enclosed mass derived from our models. 

Pastoriza et al. (1999), using ground-based gas kinematics, 
assume that the gas disk in NGC 3379 is nearly face-on. Their 
mass estimate inside 1.3" is 7 x 1O 9 M , strongly contradicting 
our result. However, their large mass is due mainly to their disk 
inclination assumption. If one uses the inclination as measured 
from the dust ring, then the enclosed mass closely equals our 
value of 1 x 10 8 M Q ; only with HST can we precisely determine 
the gaseous disk configuration and determine whether the gas is 
even in a disk at all. The HST/FOS spectrum includes the [Nil] 
emissiom, although since it is a single pointing, we do not have 
the radial information to extract a meaningful dynamical anal- 
ysis from the gas. However, the width of the line does provide 
some information on the underlying mass. The measured a is 
200km s" 1 , in agreement with the mass deduced here; using the 
mass from Pastoriza et al. would imply a a around 700km s" 1 . 
To further understand the gas kinematics, we must have spatial 
information. 

Van der Marel et al. (1990) apply 2-integral flattened models 
for NGC 3379 using constant anisotropy. They find an accept- 
able fit with an inclination of 60° and anisotropic orbits. They 
suffer from the same problem as we do for the inclination esti- 
mate; without inclusion of a dark halo it is not possible to obtain 
adequate inclination constraints. Furthermore, since our ax- 
isymmetric models are fully general while their models specify 
a form for the distribution function, we cannot directly compare 
our results to theirs. Van der Marel et al. conclude that there is 
no evidence to support the need for a third integral using data 
from the first and second moments of the LOSVD from Davies 
& Birkinshaw (1988). Results presented here, which are based 
on higher S/N data and are determined from the full LOSVD, 
show that the best model requires three integrals and is not con- 
sistent with a 2-integral distribution function; the internal dis- 
persions along the long axis in the radial and 9 directions are 
significantly different, in direct conflict with a 2-integral distri- 
bution function. As seen in Fig. 13, the ratio of the 9 to the 
radial dispersion in our best model have a radial variation from 
0.6 to 1.2, with an average around 0.75. For comparison, this 
ratio is around 0.6 for the spheroid component of our Galaxy 
using RR Lyrae halo stars (Layden 1995). 

In addition to the measured black-hole mass, the orbital 
structure appears to be a robust feature as well. The radial 
anisotropy profile along the major axis in Fig. 13 shows that 
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this profile is nearly independent of input inclination. For all ac- 
ceptable models, the velocity ellipsoid is tangentially biased in 
the central bin and radially biased in the mid-range radii. These 
orbital properties are generally true along any position angle. 
This central tangential bias appears in other 3-integral models 
studied (Gebhardt et al. 2000) and may be a common feature 
of early-type galaxies. Tangential bias in the central regions 
around black holes appear to be a common feature in simula- 
tions (Quinlan et al. 1995, Quinlan & Hernquist 1997, Merritt 
& Quinlan (1998), and Nakano & Makino 1999). The radial 
bias in the mid-range radii has also been seen before; using 3- 
integral models, both Gerhard et al. (1999) for NGC 1600 and 
Saglia et al. (1999) for NGC 1399 find radial motion through- 
out the main body of the galaxy. Thus, we may be begin- 
ning to measure common features of the orbital properties in 
ellipticals — tangential motion in the central regions and radial 
motion at mid-range radii. At large radii, we must include dark 
halos since they make a significant contribution there. The next 
step is to determine dominant evolutionary effects from theoret- 
ical models and will require detailed comparisons with N-body 
simulations. 

Magorrian et al. (1998), using the same ground-based data, 
fit two-integral models to the observed second-order velocity 
moment profiles and require a BH of 3. 9+Q4 x 10 8 M q . Our 
more general three-integral models would allow an even larger 
range of BH masses had we similarly used only the LOS VD's 
second moments. When we fit to the full LOSVDs, how- 
ever, their shape requires that our models be mildly radially 
anisotropic, reducing the required BH mass by a factor of at 
least two over Magorrian et aVs result. The models presented 
in this paper still contain restrictive assumptions — namely con- 
stant M/Ly and spheroidal symmetry — requiring yet more mod- 
els to study black hole properties in complete generality. 

Kormendy (1993) and Kormendy & Richstone (1995) find 
that there is a correlation between the black hole mass and the 
mass of the spheroid. Whether this relation is a real correlation 
or an upper envelope that extends to smaller black-hole masses 
remains to be seen. As the data quality and the modeling tech- 
niques improve, we should be better able to constrain the re- 
lationship between black hole and spheroid mass, and begin to 
measure correlations with the orbital characteristics. 
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NGC 3379 HST Photometry Results 
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